Hydrodynamical model for J/ip suppression and elliptic flow 



A. K. ChaudhurE 

Variable Energy Cyclotron Centre, 
1/AF, Bidhan Nagar, Kolkata 700 064, India 

In a hydrodynamic model, we have studied J/ip suppression and elliptic flow in Au+Au collisions 
at RHIC energy y / s=200 GeV. At the initial time, J/ip's are randomly distributed in the fluid. 
As the fluid evolve in time, the free streaming J/ip's are dissolved if the local fluid temperature 
exceeds a melting temperature Tj/^. Sequential melting of charmonium states (xc, ipl and J/ip), 
with melting temperatures T Xc = T^,, « 1.2T C , Tj/,/, ~ 2T C and feed-down fraction F ~ 0.3, is 
consistent with the PHENIX data on J/ip suppression and near zero elliptic flow for J ftp's. It is 
also shown that the model will require substantial regeneration of charmoniums, if the charmonium 
states dissolve at temperature close to the critical temperature, T Xc — T^,, < T c , Tj/^ « 1.2T C . The 
regenerated charmoniums will have positive elliptic flow. 

PACS numbers: PACS numbers: 25.75.-q, 25.75.Dw 



Matsui and Satz |l|, predicted that in presence of 
quark-gluon plasma (QGP), due to color screening, bind- 
ing of a cc pair into a J/ip meson will be hindered, leading 
to the so called J/ip suppression in heavy ion collisions. 
Over the years, several groups have measured the J/ip 
yield in heavy ion collisions (for a review of the data 
prior to RHIC energy collisions, and the interpretations 
see Refs. [1, 9). In brief, experimental data do show 
suppression. However, this could be attributed to the 
conventional nuclear absorption, also present in pA colli- 
sions. At RHIC energy (^1=200 GeV), PHENIX collab- 
oration has made systematic measurements of J/ip pro- 
duction in nuclear collisions. They have measured J/ip 
yield in p+p collisions at RHIC and obtained the refer- 
ence for the basic invariant yield 0, [H, 0. Measurements 
of J/ip production in d+Au collisions H, 0] give refer- 
ence for cold nuclear matter effects. J/ip production in 
d+Au collisions are consistent with cold nuclear matter 
effect quantified in a Glauber model of nuclear absorp- 
tion with a a bs — 2 ± 1 mb (8j. Cold and hot nuclear 
matter effects are studied in Au+Au and Cu+Cu colli- 
sions, where yields are measured as a function of collision 
centrality At RHIC energy, it has been argued 

that rather than suppression, charmonium's will be en- 
hanced [H, EH ■ Due to large initial energy, large number 
of cc pairs will be produced in initial hard scatterings. 
Recombination of cc can occur enhancing the charmo- 
nium production. Apparently, both the PHENIX data on 
J/ip production in Au+Au and in Cu+Cu collisions, are 
not consistent with models which predict J/ip enhance- 
ment [Hjdlj]. J/ip's are suppressed both in Au+Au and 
Cu+Cu collisions, suppression is more in central colli- 
sions than in peripheral collisions. 

Recently PHENIX collaboration [lj[ measured the el- 
liptic flow for J/ip in 20-60% Au+Au collisions. Ellip- 
tic flow for J/ip is an important observable. It can test 
whether or not the J/ip production is dominated by re- 
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combination. In p+p and Au+Au collisions, PHENIX 
has measured semi-leptonic decay electrons from heavy 
quarks [15]. In Au+Au collisions decay electrons has pos- 
itive elliptic flow. If recombination of cc is a major source 
of J/ip in Au+Au collisions, as suggested in [l2|, [r|, 
J/ip's will inherit some of their flow. The PHENIX mea- 
surements of J/ip elliptic flow has large error bars. Inte- 
grated V2 is consistent with zero, v 2 = —0.10+0.10+0.02. 

Recently, in a hydrodynamic based model [l(| [l?], OH , 
centrality dependence of J/ip suppression in Au+Au and 
Cu+Cu collisions are studied. In the model, J/ip's are 
randomly produced in initial NN collisions. As the 
fluid evolve, free streaming J/ip's are melted if the lo- 
cal fluid temperature exceed a threshold temperature 
Tj/^. Sequential melting of charmonium states (\c, ip 1 
and J/ip), with melting temperatures T Xc = T^„ w 1.2T C , 
Tj/^ ~ 2T C and feed-down fraction F ss 0.3, explains the 
PHENIX data on the centrality dependence of J/ip sup- 
pression in Au+Au collisions. J/ip px spectra and the 
nuclear modification factor in Au+Au collisions are also 
well explained in the model. The model leaves little or 
no room for regeneration of J/ip due to recombination of 
cc pairs, as suggested in [l2|, [r| . The model is also con- 
sistent with zero elliptic flow for J/ip's. Initially J/ip's 
are produced randomly, they do not have any flow. In 
later times also, the free streaming J/ip's can not acquire 
any flow. 

The dissociation or melting temperatures of different 
charmonium states obtained in (l6l . UtI HH are in agree- 
ment with potential model calculations for charmonium 
states at finite temperature 19]. Due to heavy mass of 
charm quarks, charmonium states can be studied in non- 
relativistic potential models. Ground state properties of 
charmonium stated are well explained using the Cornell 
potential, v(r) = ar — 2 with string tension a w0.2 
GeV 2 and gauge coupling a ss 7r/12 [l9|. Lattice QCD 
can provide for the heavy quark potential at finite tem- 
perature. Finite temperature potential models indicate 
that charmonium states, J/ip(lS), x c (lP) and ipr(2S) 
dissociate respectively at temperatures, Tj/^ 2T C , 
T Xc = T^i w 1 — 1.2T C 19]. However, very recently, in 
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FIG. 1: The open circles are hydro-dynamic model predictions 
for J/ip survival probability in (a) 0-5% and (b) 20-30% cen- 
trality Au+Au collisions, as a function of J/ip melting tem- 
perature (Tj/^). The shaded regions in panel (a) and (b) 
indicates PHENIX measurements for J/ip survival probabil- 
ity [iol l25l [26l] . Note that if J/ip dissociation temperature is 
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1.2T C , experimental data will be under explained. 



20, 21] quarkonia spectral function in QGP is determined 
using a potential motivated by lattice QCD results on free 
energy of static quark-antiquark pair. Surprisingly, char- 
moniums are found to melt at much lower temperature, 
Tj/t)) ~ 1-2T C , T Xc = T^i < l.2T c . 

In the present paper, we show that if the charmo- 
nium states dissolve at temperature close to the critical 
temperature, the PHENIX data on J/ip suppression in 
Au+Au collisions are not explained in the hydrodynamic 
model 

EE EEIIi- J/ip's will be more suppressed than 
in the experiment. Data on J/ip suppression can only be 
explained with substantial regeneration of J/ip during 
the evolution. It is also shown that regenerated charmo- 
niums will have positive elliptic flow. Conversely, small 
positive elliptic flow for J/ip will indicate regeneration of 
charmoniums in RHIC energy collisions. 

Details of the hydrodynamic model used here can be 
found in E3j EH- Briefly, it is assumed that in high 
energy nuclear collisions, a deconfined phase (QGP) is 
produced, which expands, cools, undergoes 1st order 
phase transition to hadronic fluid at the critical tem- 
perature (T c =164 MeV) and then further cools to freeze- 
out at temperature Tp=130 MeV. Assuming longitudinal 
boost-invariance, the space-time evolution of the fluid is 
obtained by solving the energy-momentum conservation 
equation dJT**" = , with initial conditions as deter- 
mined in [22], e.g. initial time t;=0.6 fm, initial cen- 
tral entropy density Si n i=llQ /to -3 , corresponding to 
energy density £i «30 GeV/ fm 3 . To obtain the sur- 
vival probability of J/^'s in an expanding medium, we 



proceed as follows: at the initial time r, = 0.6/to/c, 
we randomly distribute a fixed number of J /ip J s in the 
transverse plane. They are assumed to be free stream- 
ing unless dissolved in the medium. Each J/ip is char- 
acterised by 4 random numbers. Two random numbers 
(i?i,i?2) indicate its transverse position (r^), and two 
random numbers (i?3,i?4) its transverse momentum px- 
Random numbers R\ and R2 are distributed according 
to the transverse profile of the number of binary colli- 
sions (N co u), calculated in a Glauber model. Random 
numbers R3 is distributed according the power law [23| . 
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with A = 4.23 and B = 4.1, which well describe the 
invariant distribution of measured J/ip's in p+p collisions 
. The random number R4 is distributed uniformly within 
[0-2tt]. 

The survival probability of a J/ip inside the expanding 
QGP is calculated as EE]> 



Sj/^{t) = exp 



T dis {T(r ± {T')))dT' 



(2) 



where T(r±) is the temperature of the fluid at the trans- 
verse position r±, T d i S (T) is the decay width of J/ip at 
temperature T. n is the initial time for hydrodynamic 
evolution. We continue the evolution till the freeze-out 
time tf , corresponding freeze-out temperature Tf=130 
MeV. For the decay width Tdis we use, 



Fdis(T) = 00; T>Tj/^ 

r dis (T) = a(T/T c -l) 2 ; T <T JN , 



(3) 



In Eq[3l a is the thermal width of the state at T/T c =2. 
NLO perturbative calculations suggest that a > OAGeV 
[H]. Presently we use a=0A El- 

In FigQJa) and (b), hydrodynamic model predictions 
for J/ip survival probability (Sj/^,), in 0-5% and 20-30% 
centrality Au+Au collisions are shown as a function of 
the melting temperature. As expected, survival prob- 
ability increases as the melting temperature increase. 
Survival probability continue to increase till Tj/^ ex- 
ceed the peak temperature of the fluid, beyond which, 
even if Tj/^ is increases, survival probability remain un- 
changed. The shaded region in Fig[TJa) and (b) repre- 
sent the experimental survival probability as determined 
in the PHENIX experiment llfji . The cold nuclear mat- 
ter effect is subtracted out [lOl [25l [26j . It is evident 
from Fig[T] that if J/t/j's dissolve at temperatures close 
to the critical temperature , Tj /ip ~ 1-2T C , as predicted 
in recent calculations [13, [Hf, PHENIX data on J/ip 
suppression is not explained in the model. Most of the 
initially produced J/i/i's (~ 98%) are dissolved in the 
medium. The data could only be explained, if during the 
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FIG. 2: Black diamonds are hydrodynamic model [3] pre- 
dictions for J/tp elliptic flow in 20-60% Au+Au collisions. 
Elliptic flow for J ftp's, randomly generated according to the 
thermal distribution in EqU in 0-5% and 20-30% and 20-60% 
Au+Au collisions are shown as blank circles, squares and tri- 
angles. The black circles are PHENIX measurements [l4[ for 
J/tp elliptic flow in 20-60% Au+Au collisions . 



evolution, cc pairs recombine to regenerate J/tp. In the 
discussion, we have neglected the effect of higher states 
Xc and tpi. Nearly 30-40% of observed J ftp's are from 
decay of the higher states. However, since Xc and tp/ 
melt at temperature lower than that for J/tp's, their in- 
clusion will require further regeneration of charmoniums. 
If on the other hand, J/tp's do survive high tempera- 
ture (e.g. Tj/tp ~ 2T C ), experimental data is explained 
in the model, without any need of recombination of cc 
pairs. Indeed, as shown in 18], PHENIX data on cen- 
trality dependence of J ftp suppression are well explained 
in the model with sequential melting of charmoniums 
states with melting temperatures, T Xe — T^i w 1.2T C , 
Tj/ip ~ 2T C and feed-down fraction F » 0.3, with little 
or no scope for recombination of cc pairs and regeneration 
of charmoniums in the QGP phase. 

Hydrodynamic model of J ftp suppression thus indi- 
cate that if the charmonium states dissolve at temper- 
ature close to critical temperature, more than ~ 98% 
of observed J/ip are regenerated from cc recombination 
during the QGP phase. If on the other hand, charmo- 
nium ground state survive high temperature, regener- 
ation of charmoniums will be minimum. Regenerated 
charmoniums can have positive elliptic flow. As men- 
tioned earlier, PHENIX collaboration measured single 
electrons from the semi-leptonic decay of heavy flavors 
in p+p and Au+Au collisions [l5| . Semi-leptonic decay 
electrons in minimum bias Au+Au collisions have posi- 
tive elliptic flow, indicating positive elliptic flow for the 
parent heavy quarks. If J/tp's are produced from cc re- 
combination in the QGP phase, they could inherit their 
flow. To obtain an idea about the elliptic flow for J/ip, 
due to recombination, we randomly generate J/tp's at 



the critical temperature T c according to the equilibrium 
distribution, 



E d?p ^ ex P(~P- u ( x )/ T c) ( 4 ) 

where u(x) is fluid velocity at the freeze-out surface at 
T = T c . Eq[4] is certainly an assumption. J/tp's can be 
regenerated throughout the deconfined phase. However, 
EqH] suffice to give the general idea of flow for the regen- 
erated J/tp's. Fluid velocity distribution at the critical 
temperature T c can be obtained from the hydrodynami- 
cal model of evolution. In Fig|2l differential elliptic flow 
for J/tp's, randomly generated according to the distribu- 
tion in EqH in 0-5% and 20-30% and 20-60% Au+Au 
collisions are shown. While in central (0-5%) collisions, 
elliptic flow is small, in peripheral (20-30% or in 20-60%) 
Au+Au collisions, regenerated charmoniums have sub- 
stantial elliptic flow. For example J/tp's produced ran- 
domly according to the distribution Eq|4]have ~ 15-20% 
v 2 at p T =2 GeV. In Figd black circles are PHENIX mea- 
surements for J ftp elliptic flow in 20-60% Au+Au colli- 
sions. The error bars are large and definitive conclusion 
about J/tp elliptic flow can not be drawn. Apparently, 
regenerated charmions have more elliptic flow than in ex- 
periment. In Fig[2l we have also shown the elliptic flow 
for J/tp in the hydrodynamic model (l8| , with melting 
temperatures, Tji^ = 2T C , T Xc = T^, = l.2T c and feed 
down fraction F=0.3 (the black circles). As stated ear- 
lier, the model is consistent with PHENIX data on J/tp 
suppression without any recombination. The elliptic flow 
in the hydrodynamic model is consistent with zero flow. 
It is expected also. In the model, J/tp's are randomly 
generated within the angular range (0 — 2ir). Initially 
they have zero elliptic flow. During the evolution also, 
free streaming J/tp's do not acquire flow. 

Before we summarise, we note that it is possible that 
difference between elliptic flow of regenerated J/tp's and 
J/tp's surviving the hydrodynamic evolution, may not be 
as large as depicted in FigO We have assumed J/tp's are 
free streaming. J/tp's are formed early in the collision. 
The formation process may induce elliptic flow. Then el- 
liptic flow of J/tp's surviving hydrodynamical evolution 
will be larger than obtained presently with the assump- 
tion of free steaming. It is difficult to ascertain the in- 
duced elliptic flow in initial J/tp's and uncertainty at the 
present evaluation can not be quantified. Present model, 
with the assumption of free streaming give the lower limit 
of elliptic flow of J/tp's surviving the hydrodynamic evo- 
lution. Also it is unlikely that regenerated J/tp will be 
equilibrated at T c . J/tp's are massive. Their mass is 
much larger than T c . Elliptic flow of non-equilibrated or 
partially equilibrated J/tp will be less than that for com- 
pletely equilibrated J/tp's. Present evaluation, with the 
assumption of equilibration then give the upper limit of 
elliptic flow of regenerated J/tp's. 

To summarise, in a hydrodynamical model, we have 
studied J/tp suppression in Au+Au collisions at RHIC. 
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In the model, at the initial time, J/tp's are randomly 
distributed in the fluid. As the fluid evolve in time, 
the free streaming J ftp's are dissolved if the local fluid 
temperature exceeds a threshold temperature Tj uu. It is 
shown that if, as suggested in some recent works |20l l2l| , 
J ftp do not survive much above the critical temperature, 
PHENIX data on J ftp suppression are not explained in 
the model. Most of the initially produced J/tp's are dis- 
solved in the medium and the model predictions severely 
underestimate the PHENIX data on J ftp suppression in 
Au+Au collisions. The data will demand recombination 



of cc pairs in the QGP phase. In a simple model, we have 
also shown that the regenerated J ftp will have positive 
elliptic flow. On the other hand, if J/tp's can survive at 
high temperature, (e.g. Tj/^> ~ 2T C ), sequential melt- 
ing of charmonium states (%c, tpf and J/tp), with melting 
temperatures T Xc = T^, w 1.2T C , Tj/^ w 2T C and feed- 
down fraction F w 0.3, is consistent with PHENIX data 
on J/tp suppression. There is no scope for recombination 
of cc pairs. The model also predict zero elliptic flow for 
J/tp's. 
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